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ABSTRACT 

Gravitational scattering between massive planets has been invoked to explain the eccentricity dis- 
tribution of extrasolar planets. For scattering to occur, the planets must either form in - or migrate 
into - an unstable configuration. In either case, it is likely that a residual gas disk, with a mass 
comparable to that of the planets, will be present when scattering occurs. Using explicit hydrody- 
namic simulations, we study the impact of gas disks on the outcome of two-planet scattering. We 
assume a specific model in which the planets are driven toward instability by gravitational torques 
from an outer low mass disk. We find that the accretion of mass and angular momentum that occurs 
when a scattered planet impacts the disk can strongly influence the subsequent dynamics by reduc- 
ing the number of close encounters. The eccentricity of the innermost surviving planet at the epoch 
when the system becomes Hill stable is not substantially altered from the gas-free case, but the outer 
planet is circularized by its interaction with the disk. The signature of scattering initiated by gas 
disk migration is thus a high fraction of low eccentricity planets at larger radii accompanying known 
eccentric planets. Subsequent secular evolution of the two planets in the presence of damping can 
further damp both eccentricities, and tends to push systems away from apsidal alignment and toward 
ant i- alignment. We note that the late burst of accretion when the outer planet impacts the disk is 
in principle observable, probably via detection of a strong near-IR excess in systems with otherwise 
weak disk and stellar accretion signatures. 

Subject headings: accretion, accretion disks — scattering — planets and satellites: formation - 
planetary systems: formation — planetary systems: protoplanetary disks 



1. MOTIVATION 

The eccentricity distribution of massive extrasolar 
planets beyond the tidal circularization ra dius is broad, 
with a median eccentricity < e >~ 0.25 (jMarcv et alJ 
20050 . Since giant planet formation models predict 
that the envelope is accreted when the planet has 
an almost circular orbit, the observed eccentricities 
are generally interpreted as evidence of evolutionary 
processes that occur subsequently. One possibility 
is that the same gravitational torques from the gas 
disk t hat drive orbital migration also excite eccen- 
tricity (lArtvmow cz et~aTlll991t lOgilvie fc Lubowl l2003t 



iGoldreich fc Saril I2003t iMoorhead fe Adamsl 120081) . Al- 
ternatively, the eccentricity may arise from gravita- 
tional instabil i ty wi t hin an unstable multiple system 
dRasio fe Fordl 119961: IWeidenschilling fc Marzarl ' Il996t 
iLin fc Idalll997l ). A-body experiments show that a wide 
variety of unstable initial conditions - which can in- 
clude systems with two, three or many giant planets - 
evolve under scattering to yield an eccentricity distri- 
bution similar to that o bserved ([Chatter jee et al.l 120071 : 
Uuric fc Trema ine 2007). The agreement between the- 
ory and observations supports the scattering scenario, 
though other possibilities are not ruled out. 

To date, most studies of planet-planet scattering have 
considered purely gravitational evolution from gas-free 
initial conditions. Gas has been ignored primarily for 
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reasons of computational necessity, since fairly robust 
arguments suggest that at least small amounts of gas 
are likely to survive to the epoch when scattering occurs 
(planetesimal disks will also be present, though these will 
likely be dynamically important only in systems similar 
to the Solar System, in which at least some giant planets 
at large radii have low envelope masses). In particular, 
for two-planet systems most initial separatio ns within 
the b oundary that is Hill stable for all time (jGladmanl 
are unstable on very short time scales. Realisti- 
cally, an unstable two planet system is therefore likely to 
arise only if the planets are driven toward such a state, 
for example due to gas disk migration. Efficient migra- 
tion requires a disk mass that is of the same order as the 
planet mass. Multiple planet systems with N p > 2 ex- 
hibit broader instability regions, but even in this case the 
most likely scenario is that instability sets in once eccen- 
tricity damping due to the presence of gas becomes small 
enough. Again, this will occur once (to order of magni- 
tude) the gas disk mass falls to the level of the planet 
mass, and scattering is likely to occur in an environment 
that is gas poor but not gas free. 

In this paper we use hydrodynamic simulations to 
study the influence of low mass residual gas disks on 
the outcome of gravitational scattering. As noted above, 
there are many possible scenarios that yield different 
initial conditions for scattering experiments - here we 
adopt a simple model that is computationally tractable. 
We assume that two Jupiter-mass planets form on well- 
separated circular orbits. Models of giant planet forma- 
tion suggest that the onset of runaway accretion occurs 
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TABLE 1 

Simulation parameters 



Case 


(AU) 


m out j 
(Mj) 


(10 32 erg s" 1 ) 


< Lace > 

(10 32 erg s' 1 ) 


C/ 

{Ccrit ) 


A 


0.00094 


1.48 


6.26 


1.76 


1.027 


B 


0.0075 


1.54 


11.37 


2.54 


1.040 


C 


0.044 


1.40 


11.51 


2.75 


1.005 


D 


0.11 


1.52 


18.82 


4.68 


1.034 


E 


0.30 


1.61 


18.36 


5.82 


1.045 



Note. — The closest approach during the initial scattering is 
d-close- Subscript / refers to the quantity at the end of the simu- 
lations, 10 4 yr. < Lace > is the mean accretion luminosity for the 
first 500 years of the simulation. 



when the cores attain a mass of 10 M e (Hubickv i et al.1 
|2005| ). and, since neighboring cores are lik ely to form with 
a sep aration of ~ 10 mutual Hill radii (jKokubo fe Ida! 
1998), this suggests that plausible initial separation ra- 
tios are ~ 1.3 or greater. If the gas disk interior to 
the outer planet is relatively weak migration driven by 
the residual outer gas results in convergence of the or- 
bits. The initial separation alluded to above lies close 
to the 3:2 resonance, so depending upon the precise ini- 
tial separation convergent migration might result in res- 
onance crossing. If the resonance is encoun tered, the re- 
sult can either be trapping into resonan ce jLee fc Peald 
120021 : iSnellgrove et_alH2001fc iKlev etTaTlteOOSl ) 4 or. if the 
disk is turbulent ([Adams et al. 2008) , avoidance of trapp- 
ping and ongoing convergence of the orbits. In the latter 
case the outer planet continues to migrate inward on a 
viscous time scale until the separation of the planets is 
driven past the Hill stability criterion, and a scattering 
event occurs. The goal of this paper is to answer the 
question: what effect does the remnant disk now have on 
the scattering planetary system? 

2. SIMULATION METHOD AND INITIAL CONDITIONS 

Hydrodynamic simulations of the disk over dynami- 
cally interesting time scales are computationally difficult. 
In this preliminary study we consider only a small num- 
ber of cases, and set up the initial conditions so that 
scattering occurs early in the simulation. Using purely 
A^-body runs, we first identified a large number of initial 
conditions in which two Jupiter mass planets, orbiting 
a Solar mass star on nearly circular orbits separated by 
2.5 to 3 mutual Hill radii, are about to experie nce a first 
close encounter. For realistic migration rates (jArmitagd 
l2007f ) instability is likely to set in at approximately this 
separation. The orbital radius of the inner planet was set 
to 3 AU. The planets were nearly coplanar, with small 
inclinations i < 1°. We chose five initial conditions from 
these Af-body runs spanning a range in the closest ap- 
proach distance attained during the first encounter, sum- 
marized in Table [1] 

We evolve the initial conditions for 10 4 yr both as 
a purely A^-body system, and hydrodynamically includ- 
ing the influence of an external gas disk. The rem- 
nant disk is assumed to extend from 5-10 AU with a 
surface density profile £(r) = Eo^ ~ 3 ^ 2 , normalized so 
that the gas disk contains one Jupiter mass of gas. 
This setup is quite similar to the scenario considered 

4 The existence of a number of known resonant pairs of planets 
supports the idea that trapping occurs in at least some cases. 
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Fig. 1. — An example snapshot from case C illustrating the 
interaction between t he outer pla net and the disk. Plotted (using 
the program splash, Price 2007) is the column density, with the 
star (dotted circle) at the center and the two planets (filled circles) 
in orbit. The outer planet has induced density waves in the disk, 
from which it is accreting gas. This snapshot is ~ 30 years after 
the close encounter between the planets. 



bv lMoorhead fc Adamsl (|2005l ). in which two-planet sys- 
tems that had cleared an inner gap were driven toward 
instability via Type II migration of the outer planet. 
These authors include the effect of a remnant disk on 
planet-planet scatte ring via an eccent r icity d amping pre- 
scription. Likewise, IChatteriee et al.l (|2007l ) studied the 
damping effect of a disk on three-planet scattering by 
coupling a planet-disk torque prescription with a 1-D 
disk model. In most of our models we assume that the 
inner disk is weak enough as to be neglected. In the 
regime we consider considerable depletion of the inner 
disk is likely, since the angular momentum of the planets 
is large enough to significantly retard the inflow of the 
outer disk while the inner disk continues to accrete on 
the normal viscous time scale. However, full depletion 
of the inner disk is only an approximation 5 , and recent 
calculations suggest that some gas will either survive in 
the inner region or flow past t he planets from the outer 
disk to resupply the inner disk (jCrida fc Morbidellill20"07l ; 
iLubow fc D'Angelol [20061 ). We therefore ran one addi- 
tional simulation to ascertain the effect of a small mass 
of gas int erior to the planets on the results. Using the 
results of ILubow fc D'Angelol (|2006f ), we estimated that 
reasonable amounts of material in the depleted inner disk 
will likely be of the order 10% of the initial disk profile 
([Lubow fc D'Angelol 120061 ). so that the inner disk has 
E(r) = 0.1£ r- 3/2 . If the inner disk extends from 0.75 - 
2 AU, it contains ~ 0.06M,/ of material, much less than 
the outer disk. If the inner planet scatters into this in- 
ner remnant a small amount of eccentricity damping and 
mass loading onto the planet could occur. 

We evolve the system using the GAD GET-2 smoothe d 
particle hydrodynamics (SPH) code (Springcl 2005), 
modified so that the star and planets behave as sink 
particles (Bat e et al.lll995l) . The accretion radius of the 



5 Unless ma gnetic effects Chiang fc Murray-Clay 2007) or pho- 
toevaporation IjClarke et alj 20011 ) operate and clear a true cavity 
out as far as the radius of the inner planet. 
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Fig. 2. — The minimum separation (equation [T} for each case 
in units of mutual Hill radii. Thick: the simulations with the gas 
disk. Thin: the same initiaj conditions run with no disk. • , , 

planets is set to ~ ouRj. rne star and planets interact 
with each other and the gas particles only via gravity, 
which is effectively unsoftened. A further modification 
to the code was made to ensure that the force between 
the star and planets is computed directly, and not with 
a tree-code approximation. For the runs presented here 
the disk is represented by ~ 1.75 x 10 5 particles. At this 
resolution the Hill radius of the planets is resolved by 
several SPH smoothing lengths throughout the disk. We 
performed a resolution test by placing a single planet on 
an eccentric orbit that initially made excursions into the 
disk at apastron. The evolution of the planet's mass, ec- 
centricity, and semi-major axis showed good agreement 
between resolutions of 1.5 x 10 5 and 3.0 x 10 5 particles. 
Figure [1] shows an example of one of our runs, showing 
the interaction of the outer planet with the disk approx- 
imately 30 years after the scattering. 

3. RESULTS 

Following the scattering event, we evolve the system 
hydrodynamically up to the time when the separation of 
the planets is great enough that - in the absence of gas 
- they would be Hill stable. This requires a relatively 
short integration of the order of 10 4 yr in duration. It 
would be desirable to continue the hydrodynamic sim- 
ulation until all of the gas has been dispersed, but this 
is both computationally more difficult and would require 
the addition of additional physics (such as photoevapo- 
ration) describing the final clean-up of the gas disk. In 
the paper, we instead study the longer term evolution 
using N-body calculations that include an approximate 
treatment of the effects of gas disk damping on planetary 
eccentricity. 

3.1. Short-term evolution 

Figure [5] shows how the separation of the planets 
evolves with time in runs with and without gas. As a 
measure of the stability of the systems, we calculate the 
minimum separation, defined by the instantaneous val- 
ues of the semi-major axes a and eccentricities e of the 
inner and outer planets by 



a / AU 

Fig. 3. — The mean values of the eccentricities and semi-major 
axes for all the simulations. Filled symbols: simulations with disks. 
Open symbols: simulations with no disks. Arrows connect the same 
initial conditions with and without disks for the outer planets. 

Figure [2] shows this quantity as a fraction of the mutual 
Hill radius 

where a is the average of the planet's semi-major axes 
and M* is the stellar mass. For nearly-circular planets, 
the system is Hill sta ble when the o rbital parameters 
satisfy A min > 3.5Rh (|Gladmanlll993f) . Though the or- 
bits in these simulations arc moderately eccentric, we use 
this measure as a intuit ive guide to short-term stability. 
iMarchal fc Bozisl (|1982f ) give a general condition for Hill 
stability, parameterized by a combination of the bodies' 
masses and the system's angular momentum and energy, 
£ oc L 2 E. A critical value ( cr u can be found, above which 
the system is stable. In the non-dissipative case L and E 
are integrals of motion, and all of our simulations begin 
below the critical value. With dissipation in the form of 
the gas disk, the energy and angular momentum of the 
three bodies are no longer constant, and £ can evolve. In 
table Q] we give the value of ( as a fraction of the critical 
value at the end of the simulations. All of the dissipative 
cases have achieved Hill stability at 10 4 yr. 

When no gas is present, the systems do not move to- 
ward stability over the timescale simulated, and all but 
one of them (case E, which is likely unstable on a some- 
what longer timescale because of its larger close approach 
distance) continue to undergo strong scattering interac- 
tions. Case A with no disk appears to be moving toward 
a more stable configuration as measured by the minimum 
separation, but the measure of Hill stability £ remains 
as it must at its initial value, in the unstable regime. 
The simulations including the remnant disk all follow the 
same general trend toward an increasing level of stability 
as measured by A TOi „/i?#-. After the first 100 years of 
evolution, none of the systems with gas present undergo 
encounters with r enc < Rh, whereas the non-dissipative 
runs continue to undergo strong interactions. 

In Figure [3] we plot the mean values of a and e over 
the span of the simulations for the different cases, with 



4 




no inner disk 
with inner disk 
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time / yr 

Fig. 4. — Comparison between case E with and without an inner 
disk. Top: the minimum separation (equation [TJ for the two cases. 
Both have reached Hill stability by 5000 yr. Bottom two panels: 
the eccentricities of both planets for each case. The two planets are 
identified by thick and thin lines. The details of the eccentricity 
evolution differ when an inner disk is added, though the final result 
in the sense of stability is very similar. 

and without the gas disk. While the inner planets show 
no clear separation in distribution, there is a consistent 
difference in eccentricity for the outer planets. The ar- 
rows connect the properties of the outer planet for each 
case, with and without remnant disks. In each instance, 
the simulation with a disk ends up at a lower eccentric- 
ity. The effect of the disk on the inclination evolution of 
the systems is much less clear, with no consistent trend 
appearing in our simulations. 

We ran case E with the addition of an inner disk as 
well. Case E was chosen because its inner planet had the 
smallest periastron post-scattering, so that any effect of 
accretion from the inner disk would be maximized. The 
total angular momentum of the inner disk we simulated 
is ^4% of the inner planet's orbital angular momentum, 
and we would expect the change in eccentricity associ- 
ated with accreting this material to be small. In figure 
[4] we compare the evolution of the orbits with the dif- 
ferent disk profiles. The details of the semi-major axis 
and eccentricity evolution between the two realizations 
are different, as the initial evolution of the inner planet 
is changed and sets the system on an altered path. The 
separation measure A. m i n shows a very similar evolution, 
however, with the same trend toward stability. As in the 
case with no inner disk, the outer planet accretes signif- 
icantly from the exterior disk and circularizes. At 5000 
years after the close encounter, the system has reached 
Hill stability (C/Ccrit = 1-033) and the inner disk has 
been completely disrupted. The further evolution of the 
system will be driven by accretion onto the outer planet. 
We conclude that the presence of a depleted inner disk 
is unlikely to significantly affect the results of our simu- 
lations. 

In contrast to dynamical friction, resonant interaction 
with the gas disk, or ejection of planetesimals, the stabi- 
lizing influence of these gas disks is due to the accretion 
of material by the outer planet. As it makes excursions 
into the disk near its apastron, it is accreting material 
with higher specific angular momentum, which causes its 
eccentricity to decrease. This effect is therefore depen- 
dent on enough material being present in the disk. Tests 
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Fig. 5. — The accretion rate of the outer planet for each case, 
for the first 2000 years of the simulation. Also plotted (dashed) is 
a reference accretion rate in the absence of scattering. 

with 0.1 Mj of disk material showed very little difference 
from the runs with no gas present. By the end of the 10 4 
yr run of the simulations the mass remaining in the disk 
has been reduced to ~ 0.2— 0.4Mj, while the outer planet 
has gained mass (m out ~ 1.5Mj, see Table[T]). Due to the 
expense of the hydrodynamic simulations, we are unable 
to extend the runs over disk-driven migration timescales. 
We cannot directly guarantee that further viscous evo- 
lution of the disk will not push the planets back toward 
instability, and restart the whole cycle. However, this 
possibility is unlikely. Once the disk mass drops signif- 
icantly below the planet mass the rate of migration is 
suppressed, though the exact degree depends upon the 
disk structure (e.g. Armitage 2007). The damage done 
to the outer disk by the initial scattering and subsequent 
accretion will have reduced the migration rate by a fac- 
tor of several, making it unlikely that the disk will be 
able to push the system back toward instability during 
its remaining lifetime. A second round of instability is 
even less likely given the fact that photoevaporation of a 
sub-Jupiter mass gas disk b y direct stellar illumina tion 
occurs on a rapid time scale ((Alexander et al.|[2006l ). 

The accretion levels during the outer planet's initial in- 
teractions with the disk, when its eccentricity is damped 
most rapidly, are quite high. In figure [5] we show these 
rates for the first 2000 years of the simulations, i.e. im- 
mediately after the first scattering event. The cases with 
the weakest initial encounters exhibited deeper penetra- 
tion into the disk, and have higher accretion rates than 
the cases with stronger interactions. With the exception 
of cases A and B, which show secondary bursts of accre- 
tion at ~ 1000 years, the main accretion onto the planets 
occurs over the first 500 years of the simulations. As a 
reference accretion rate, we ran a simulation with a sin- 
gle planet orbiting at the 2:1 resonance with the inner 
edge of the disk for 1000 yr. This rate is plotted in figure 
[5] as well, and has a mean value ~ 5 x 10~ 8 M Q yr . 

3.2. Long-term evolution 

We terminated our SPH simulations with the inner 
planet on an eccentric orbit (typically e ~ 0.2 — 0.3) 
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Time (yr) 

Fig. 6. — Long-term orbital evolution of two Jupiter-mass planets 
with a drag force being applied to the outer planet. Top panel: 
Eccentricities of the inner (grey line) and outer (black line) planets 
in time. Bottom panel: the relative apsidal angle, i.e. the angle 
between the longitudes of pericenter of the two planets, zo\ — vji- 



and the outer planet still interacting with the remnant 
gaseous disk, though much more weakly than immedi- 
ately after the scattering encounter. The lifetime of 
the remnant gaseous disk is uncertain: 10 5-6 yr is a 
reasonable estimat e based on models o f dissipation via 
photo evaporation (jClarke et al.l 120011 : lAlexander et al. 
2006 ) and observational constraints (|Simon &: Prate 



19951 IWolk fc Walter|[l996l ). Given that the outer planet 



is experiencing orbital damping and is also gravitation- 
ally interacting with the inner planet, an interesting 
question is whether the inner planet's eccentricity can 
be decreased. We therefore investigate the long-term dy- 
namical evolution of these systems. 

We mo del the system using the hybrid integrator 
Mercury ( Chambers! Il999f ) . The physical scenario that 
we have in mind is that the eccentricity of the outer 
planet will be damped on long time scales via resonant in- 
teraction with the remaining gas disk, once it has settled 
down following the perturbation induced by the scatter- 
ing. The exact form of the damping is not likely to be 
critical, so for simplicity we use a method that we have 
tested for other applications and apply a drag force to 
the outer planet in the form of a headwind which is pro- 
portional to the difference in velocity between the planet 
and the gas (whi ch is assumed to be p artially pressure 
supported, as in iThommes et af] l2003f ). This force is 
there fore equivalent to aerodynamic drag (|Adachi et al.l 
197(f). For our purposes we treat the strength of the 
damping as a free parameter, and calibrate the resulting 
eccentricity damping to match the damping of the outer 
planet in our SPH simulations. The effect of this force is 
to damp e on a timescale of a few times 10 years with 
little damping of a, as long as e is less than ~ 0.2. We 
therefore believe that our implementation is a reasonable 
approximation. 

Figure [6] shows the evolution of two Jupiter-mass plan- 
ets for a case that started with the inner planet with a 
= 2.5 AU and e = 0.2 and the outer planet at 5 AU 
with e = 0.1. The outer planet's eccentricity is quickly 
excited by the inner planet and the two planets undergo 
regular secular oscillations with a period of about 15,000 
years. The initial amplitude of eccentricity oscillations 
is close to the inner planet's starting eccentricity, and 
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Fig. 7. — Long-term orbital evolution of two Jupiter-mass planets 
with a time-dependent drag force applied to the outer planet. The 
drag force decreases linearly to zero after 10 years and no drag is 
applied for the rest of the simulation. 

the two values remain out of phase. The damping de- 
creases the eccentricities of both planets in time, but 
does not affect the secular frequency. At the end of the 
simulation, the inner/outer planets have median eccen- 
tricities of 0.036/0.032. The damping also affected the 
planets' apsidal aligment, as measured by the difference 
in the longitudes of pericenter, i.e., w\ — W2- The two 
planets quickly became locked in libration about apsidal 
alignment (zu% — txji = 0), but the libration amplitude 
increased in time until the nodes began to occasionally 
circulate after ~ 6 x 10 5 years. For all the configura- 
tions that we tested, evolution moved the planets away 
from apsidal alignment and toward ant i- alignment. For 
cases in which the outer planet started with e = 0, 
the system quickly tended toward anti-alignment and 
the amplitude of libration decreased steadily to from 
> 9 0deg at early times to l ess than 50deg after 1 Myr 
(sec IChiang fc Mu7ravl[200l . 

The finite lifetime of the gaseous disk may limit the 
damping of the eccentricity of the outer, and therefore 
also inner, planet. The case from figure [6] assumed con- 
stant damping from the disk over 1 Myr. However, the 
outer disk is being depleted by both accretion onto the 
outer planet and photoevaporation. Figure [7] shows the 
evolution of the same two planets from figure [6] but 
assuming the damping dissipated linearly from its de- 
fault value to zero after 10 5 years, a more reasonable 
valu e for the length of the final stage of disk d issipation 
(e.g. IClarke et al.l 120011: lAlexander et al.ll2006l ). Eccen- 
tricities of the planets in figure [7] are far less damped 
than in figure [6j the median final eccentricities of the 
inner/outer planets are 0.127/0.12. Note that for the 
case of no damping, the median final eccentricities of the 
inner/outer planets are 0.24/0.14. Thus, even a short pe- 
riod of damping can reduce eccentricities. Note also that 
the weaker damping has slightly increased the libration 
amplitude but not enough to cause the apses to circulate. 
We regard this second setup, with a limited amount of 
damping, as the more realistic scenario. The remnant 
disk mass is reduced by the initial planet-disk interac- 
tions, and estimates of the photoevaporation timescale 
are closer to 10 5 than 10 6 yr. 

The known extra-solar multi-planet systems are pref- 
erentially located close to the separatrix between apsi- 
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dal circulation and libration (|Barnes fc Greenberj[200l . 
Near-separatrix behavior is characterized by one of the 
planets periodi cally reaching an e — state (e.g. 
iFord et aLll2005[ ). The strong damping experienced by 
the system in figure [6] has driven the system from libra- 
tion to a near-separatrix state, whereas the system from 
figure[7]is steadily librating. For cases which started with 
a zero eccentricity outer planet, the situation is reversed: 
simulations that experienced only a small amount of 
damping remained close to the separatrix but those that 
were strongly damped were driven to low-amplitude ap- 
sidal libration about anti- alignment. Thus, interactions 
between planets during this phase of a low-mass, remant 
gas disk may play an important role in establishing the 
final apsidal alignment of planetary systems. 

4. DISCUSSION 

The expense of the hydrodynamic simulations means 
that we have been able to run only a handful of real- 
izations of scattering in the presence of gas, and as a 
consequence it is unwise to draw definitive conclusions. 
Nonetheless we find, as expected, that remnant disks that 
have masses comparable to the masses of the giant plan- 
ets do influence the outcome of gravitational scattering, 
primarily via eccentricity damping that occurs due to ac- 
cretion when planets intrude upon the gas disk. If, as we 
have assumed, the remnant gas lies exterior to the orbit 
of the outermost planet, damping when the planet enters 
the disk tends to increase the planetary separation and 
stabilize the system against close encounters and possible 
ejection of one of the planets. The stabilization occurs 
rapidly enough that we were able to evolve the systems 
with gas up to the point where they would plausibly be 
Hill stable. 

The observed eccentricity distribution of extrasolar 
planets 6 matches the theoretically pr edicted one from 
simu l ations that do not incl ude gas ()Chatteriee et aLl 
l2007t Uuric fc Tremainel l2007t) . If gas is implicated in 
driving planetary systems toward instability, its effect on 
the subsequent dynamics cannot be so strong as to ruin 
this agreement. We find that this is plausible. Specifi- 
cally, models computed with an exterior disk whose mass 
equals that of the planets yield inner planet eccentricities 
that are comparable to those achieved in the absence of 
gas. Since radial velocity surveys favor the detection of 
planets at smaller orbital radii, the observed sample to- 
day is likely dominated by planets which would have been 
unaffected by gas during the scattering process. Con- 
versely, we find that the eccentricity of the outer planet 
is significantly damped by accretion of gas once it en- 
ter the disk. The dynamical dominance of accretion over 
other damping mechanisms is likely the most robust re- 
sult of this paper. 

Although we have only considered one specific scatter- 
ing scenario, it seems plausible that the generic effect of 

6 A detailed study of possib le selection biases in the Keck planet 
search (Cummins ct al. 20Q3) shows that the survey is complete 
for velocity amplitudes K > 20 ms* 1 and e < 0.6. The "observed" 
distribution is therefore a reasonable proxy for the true distribu- 
tion in the sense that neither low nor moderately high e objects 
are systematically missed. There may, however, be some upward 
bias to the derived eccentricities of detected planet s in cases where 
the m easurement signal to noise is relatively low ()Shen fc Turner! 
|200D . 



residual gas disks will be to increase the fraction of re- 
tained rather than ejected planets following scattering. If 
so, analysis of the abundance of multiple planets ought to 
provide the strongest constraints on models of this type. 
If the specific model for migration and scattering we have 
used were actually commonplace in extrasolar planetary 
systems, we would expect to find a high fraction of multi- 
ple planet systems in which the outer planet had a lower 
eccentricity than the inner planet. Observationally, it 
is known that there are hints for additiona l companions 
in ma ny systems that host known planets (|Wright et al.l 
l2007f ). but the quantitative statistical analysis necessary 
to place lim its on their masses and periods (analogous to 
the work of ICumming et aT]|2008l for single planets) has 
not, to our knowledge, been attempted. Our impression 
is that while our toy model - in which the hypothetical 
outer planet has the same mass as the detected inner one 
- could already be ruled out observationally, a more re- 
alistic scenario in which the outer planet was modestly 
less massive probably remains viable. More extensive 
theoretical work will be needed to quantify the extent to 
which residual gas disks increase the predicted fraction 
of multiple systems after scattering. 

On Myr timescales, the inner planet's eccentricity can 
be damped indirectly by the remnant gaseous disk via ex- 
citation of the outer planet and subsequent dissipation. 
The lifetime of the remnant disk is uncer t ain, but is prob- 
ably ~ 10 5 years llSimon fc Pratol [l995t I Work fc Waited 



flifflfl IClaike et all 120011 : lAlexander et all 120061 ). Given 
the large eccentricity of the inner planet and the rela- 
tive proximity of the outer planet, there is strong secular 
coupling and eccentricity oscillations are induced. The 
damping of the outer planet by the remnant gas disk 
tends to push the system away from apsidal alignment 
and toward anti- alignment; the amount of damping de- 
pends on the disk lifetime. Depending on the starting 
configuration, strong damping can produce systems that 
are n ear the apsidal separatrix (|Barnes fc Greenberd 
2006), circulating, or librating about anti-alignment. 
However, strong damping should not produce systems 
that librate about alignment. In cases with relatively 
little damping, the apsidal alignment is weakly affected, 
although eccentricities can be modestly reduced (see fig- 
ure [3). 

Finally, we note that the accretion rate onto the planet 
when the planet first enters the disk is rather high. This 
gas will, in reality, join a circumplanetary disk, and how 
quickly it is accreted by the planet will depend upon the 
evolution of that disk. However, it is clear that the lu- 
minosity of the planet could be large due to the implied 
late burst of accretion. Naively estimating the accretion 
luminosity as L acc ~ GM p M/R p , we find that for Jo- 
vian parameters the luminosity over the first 500 yr is 
of the order of 10 32 -10 33 erg s 1 . Since the accreted gas 
has angular momentum, the mass inflow onto the planet 
would be buffered by a circumplanetary disk, though 
only modestly since the viscous timesc ale of such disks 
is typ ically of the order of 10 3 yr or less (|Canup fc Wardl 
2002). Unlike the main envelope accretion event that oc- 
curs in core accretion models for giant planet formation, 
this late burst of accretion is in principle observable (as a 
near-IR excess without matching stellar accretion signa- 
tures) since the accretion luminosity of the circumstellar 
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disk at this epoch is small. If we assume that late ac- 
cretion is detectable for 10 3 yr, and occurs in 10-100% 
of young stars, then we would need to survey 10 4 to 10 5 
systems to have a chance of directly observing the last 
stage of giant planet formation. Surveys of this mag- 
nitude are well beyond current capabilities, but are not 
inconceivable in the future. 
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